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ABSTRACT 

Neutron stars and black holes often reside in binaries where the accretion rate 
onto the compact object varies by orders of magnitude. These "X-ray transients" 
are observed both in outburst (when the high accretion rate makes them X-ray 
bright) and quiescence (when the accretion rate is very low, or potentially zero). 
In a previous paper, we showed that the quiescent X-ray emission from three neu- 
tron star transients (Aql X-l, Cen X-4, and 4U 1608—522) were well represented 
by thermal emission from the neutron star's hydrogen atmosphere and that the 
emitting area was consistent with the whole surface. Previous black-body spec- 
tral fits (which are not accurate representations of the thermal spectrum) severely 
underestimated the true emitting area. In this paper, we fit hydrogen atmosphere 
models to the X-ray data for four neutron stars (the three from the previous pa- 
per, plus 4U 2129+47) and six black hole candidates (A0620-00, GS 2000+25, 
GS 1124-68, GS 2023+33, GRO J1655-40, and GRO J0422+32) with masses 
^ 3 M & . While the neutron stars are similar in their intrinsic X-ray spectra (that 
is, similar effective temperatures and emission area radii ~ 10 km), the spectra 
of two black hole candidates are significantly different, and the spectra of the 
remaining four are consistent with a very large parameter space that includes 
the neutron stars. The spectral differences between the neutron stars and black 
hole candidates favors the interpretation that the quiescent neutron star emission 
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is predominantly thermal emission from the neutron star surface. Higher quality 
data from Chandra, XMM, and ASTRO-E will yield a much better contrast. 

There are many transients which do not have clear neutron star characteristics 
(such as type I X-ray bursts or coherent pulsations) and where the mass of the 
compact object is not constrained. In these cases, it is ambiguous as to whether 
the compact object is a neutron star or black hole. Our work suggests that 
an X-ray spectral comparison in quiescence provides an additional means for 
distinguishing between neutron stars and black holes. The faint X-ray sources in 
globular clusters - thought to be either cataclysmic variables or quiescent neutron 
stars - are a class of objects which can be investigated in this manner. 



Subject headings: accretion, accretion disks — stars: neutron — stars: individual 
(Aql X-l, 4U 1608-522, Cen X-4, 4U 2129+47, GRO J0422+32, A0620-00, 
GS 1124-68, GRO J1655-40, GS 2000+25, GS 2023+33, EXO 0748-676) 



1. Introduction 

Many black holes and neutron stars are in binaries where a steady-state accretion disk 
(one that supplies matter to the compact object at the same rate as the mass is donated from 
the companion) is thermally unstable ( |Van Paradijs 199^ ; [King et al. 1996| ). This thermal 



instability results in a limit cycle - as in dwarf novae - with matter accumulating in the outer 
disk for ~months to decades until a thermal instability is reached (see Huang fc Wheeler 1989 



Mincshige fc Wheeler 19891 ), trig gering rapid accretion onto the compact object, and a 



substantial X-ray brightening (typically by xlO 3 or more, up to 10 38 erg s _1 ). Both neutron 
star (NS) and black hole candidate (BHC) systems exhibit these X-ray outbursts, separated 
by periods (~ months to decades) of relative quiescence (for recent reviews, see |Tanaka H 
|Lewin 1995| ; [lanaka fc Shibazaki 1996] |Uhen et al. 19971 ; |Oampana et al. 1998 b| ). 



Deep pointed X-ray observations of these transients in quiescence have found that BHCs 
are, on average, less luminous in quiescence than their NS counterparts (|B arret et al. 19961 ; 



Chen et al. 19971 |Narayan et al. 1997b| ; |Asai et al. 1998|) . To explain the higher quiescent 



luminosities (10 32 -10 33 erg s _1 ) of the transient NSs, Narayan et al. ( |1997b| ) suggested that 
matter continues to accrete during quiescence onto both objects, and that neutron stars re- 
radiate the in-falling matter's kinetic energy, while black holes swallow most of the accreted 
mass-energy. In this picture, the required quiescent accretion rate, M q , onto the compact 
object in BHC systems must be substantially greater than in NS systems. Current spectral 
modeling (for an advection dominated flow) of the few X-ray detected BHCs use M q of ~ 
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1/3 of the total mass transfer rate in the binary ( Narayan et al. 1997a ) whereas the more 



efficient energy release onto a neutron star requires that M q be smaller by 2-3 orders of 
magnitude. The cause of this difference in M q between the BHC and NS systems is not 
easily explained flMenou et al. 199TS| ). 



An alternative picture for the NS emission was put forward by Brown, Bildsten, & 
Rutledge ( |1998| ), who showed that NSs radiate a minimum luminosity, even if accretion 
completely ceases during quiescence. This minimum luminosity comes from energy deposited 
in the inner crust (at a depth of ~300m) during the large accretion events. The freshly 
accreted material compresses the inner crust and triggers nuclear reactions that deposit ~ 
MeV per accreted baryon there ( [Haensel fc Zdunik 199(J| ). This heats the NS core on a 10 4 -10 5 
yr timescale, until it reaches a steady-state temperature ^4x 10 7 ((M)/10 -11 M yr _1 ) a4 K 
( |Bildsten fc Brown 1997| ), where (M) is the time-averaged accretion rate in the binary. 
A core this hot makes the NS incandescent, at a luminosity L q = 1 MeV (M)/m p k 6 x 
lO 32 ((M)/lO" 11 M yr" 1 )erg s _1 , even after accretion halts (grown et al. 1998Q . The NS is 
then a thermal emitter in quiescence, much like a young NS. 

For both hypotheses, the energy source for the quiescent luminosities of the BHCs and 
NSs have different physical causes, making it meaningful to search for distinguishing spectral 
signatures. 



1.1. X-Ray Spectra of Quiescent Neutron Stars 



The first NS transient detected in quiescence was Cen X— 4 ( |Van Paradijs et al. 1987| ) 



More recently, quiescent X-ray spectral measurements have been made of Aql X— 1 ( |Ver-| 
bunt et al. 19941) and 4U 2129+47 flGarcia fc Callanan 1999|) with the ROSAT /PSPC; 



of EXO 0748-676 with Einstein IPC flGarcia fc Callanan 19991) ; and of Cen X-4 and 



4U 1608-522 with ASCA QAsai et al. 1996 b| ). The X-ray spectrum of Aql X-l (0.4- 
2.4 keV) was consistent with a blackbody (BB) spectrum, a bremsstrahlung spectrum, or 
a pure power-law spectrum ( |Verbunt et al. 1994|) . For 4U 1608—522, the spectrum (0.5- 



10.0 keV) was consistent with a BB (/cTbb ~ 0.2-0.3 keV), a thermal Raymond-Smith model 
(kT = 0.32^q;5 8 keV), or a very steep power-law (photon index 6^2)- Similar observations 
of Cen X— 4 with ASCA found its X-ray spectrum consistent with these same models, but 
with an additional power-law component (photon index pa 2.0) above 5.0 keV (recent obser- 
vations with BeppoSAX of Aql X— 1 in quiescence also revealed a power-law tail; |Campana 



|et al. 1998a| ). The origin of the observed power-law spectral components in Cen X— 4 and 
Aql X— 1 is not clear. While it has been suggested they may be due to magnetospheric 
accretion (|Campana et al. 1998b| ), spectral models of metallic NS atmospheres flRajagopal 
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& Romani 1996; Zavlin et al. 1996|) also predict hard tails. These warrant further observa- 



tional investigation. The observation of EXO 0748—676 shows it to be more luminous (by 
xlO — 50) than the other four NSs. 

In four of these five sources (the exception being EXO 0748—676), BB fits implied an 
emission area of radius ~ 1km, much smaller than a NS. This had little physical meaning 
however, as the emitted spectrum from a quiescent NS atmosphere with light elements at 
the photosphere is far from a blackbody. For a weakly-magnetic (B < 10 10 G) pure hydrogen 
or heliumQ atmosphere at effective temperatures kT e ^^ 0.5 keV the opacity is dominated by 
free- free transitions ( [Rajagopal fc Romani 199~6| ; [Zavlin et al. 1996] ). Because of the opacity's 



strong frequency dependence (oc z^~ 3 ), higher energy photons escape from deeper in the 
photosphere, where T > T cS ( [Pavlov fc Shibanov 1978| ; [Romani 1987| ; |Zampieri et al. 19"9~5j ). 
Spectral fits of the Wien tail (which is the only part of the spectrum sampled with current 
instruments) with a BB curve then overestimate T e ff and underestimate the emitting area, 
by as much as orders of magnitude flRajagopal fc Romani 19961 ; |Zavlin et al. 1996D PT. 



Rutledge et al. ( |1999| ) showed that fitting the spectra of quiescent NS transients with 
these models yielded emitting areas consistent with a 10 km radius NS. In Fig. [I], we compare 
the measured H atmosphere and blackbody spectral parameters for the quiescent NSs. The 
data for 4U 2129+47 is analysed here, while the other sources were analysed previously 
( |Rutledge et al. 1999|) . The emission area radii are larger from the H atmosphere spectra 



by a factor of a few to ten, and are consistent with the canonical radius of a NS. There is 
thus both observational evidence and theoretical motivation that thermal emission from a 
pure hydrogen photosphere contributes to - and perhaps dominates - the NS luminosity at 
photon energies at 0.1-1 keV. This makes possible the use of quiescent NS X-ray spectra as 
an astrophysical tool. 



6 The strong surface gravity will quickly (within ~ 10 s) stratify the atmosphere ( Alcock fc Illarionov 198(] : 
Romani 1987); for accretion rates ^ 2 x 10 _13 A/ Q yr _1 (corresponding to an accretion luminosity <J 2 x 



10 erg s ), metals will settle out of the photosphere faster than the accretion flow can supply them 
(Bildsten, Salpeter, & Wasserman 1992 ). As a result, the photosphere should be nearly pure hydrogen. 

7 Application of H atmosphere models to the isolated neutron stars in SNR PKS 1209—52 and Puppis 
A produced a source distance consistent with that measured through other means (assuming a 10 km NS 
radius), a lower surface temperature, and an X-ray measured column density that was consistent with that 
measured from the extended SNR (while the column density measured with an assumed BB spectrum was 
not consistent with other measurements; Zavlin et al. 1998 ; Zavlin et al. 1999 ). 
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1.2. Quiescent X-Ray Spectroscopy of NSs and BHCs 

Distinguishing between a stellar-mass black hole and a neutron star in an X-ray binary 
is a non-trivial observational problem. Although there are X-ray phenomena unique to NSs 
(for example, type I X-ray bursts and coherent pulsations), as yet, no X-ray phenomenon 
predicted to occur exclusively in BHs has been observed. In the absence of any distinctive 
NS properties, some X-ray transients are classified as BHCs if they display X-ray spectral 
and variability properties similar to those of other BHCs, such as ~30% rms frequency-band- 
limited variability accompanied by a hard spectrum, or 3-12 Hz quasi-periodic oscillations 
while the source has high X-ray intensity, although there are NS systems which display 
these properties as well (see Van der Klis 1995 for a review). While a statistical distinction 



between the quiescent luminosities of BHCs and NSs has been demonstrated (see refs. in 
the Introduction), there is overlap in the observed luminosities; thus, while the different 
average quiescent luminosities support the hypothesis that the (pre-classified) objects belong 
to distinct classes, quiescent luminosity cannot be used to distinguish bewteen a NS and a 
BHC on a case-by-case basis. A promising phenomenological distinction between NSs and 
BHCs is that, at X-ray luminosities >10 37 erg s" 1 , the 20-200 keV luminosity of BHCs is 
systematically higher than that of NSs. However, the physical origin of this difference is not 
known, and the weakness of a phenomenological distinction is its vulnerability to a single 
counter example ( |Barret et al. 1996| ). 



By far the most solid technique is measuring or constraining the mass of the compact 
object via radial velocity measurements of the optical companion. If the resulting optical 
mass function indicates that the compact object mass exceeds 3M , then it is likely to be 
a black hole, as the mass-limit of a NS has been calculated to be below 3.2M Q ( [Rhoades 1& 



Ruflini 1974| j |Chitre fc Hartle 1976| ). The large amount of progress in this method has given 



us several very secure black holes ( |McClintock 1998 ). 



A new spectroscopic distinction between transient NSs and BHs - based on the presence 
of a neutron star's photosphere - would be a valuable classification tool. We present the first 
such comparisons here, where we report a spectral analysis that uses an accurate emergent 
spectrum from a NS to fit the quiescent X-ray spectra of six transient BHCs with measured 
mass-functions (GRO J0422+32, A0620-00, GS 1124-68, GRO J1655-40, GS 2000+25, 
and GS 2023+33) and four transient neutron stars. 

We begin in §Q by describing the BHCs and NS's we have chosen for this comparison, 
as well as the hydrogen atmosphere models. Three of the NS are from our previous study 
( IRutledge et al. 1999Q , Aql X-l, Cen X-4, 4U 1608-522; and we discuss a fourth here, 



4U 2129+47. We show in § || that although the neutron stars occupy a narrow range of 
effective temperatures and emitting area radii, no such relation is found among the BHCs. 
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This implies that the H atmosphere spectrum may be used as a tool to distinguish between 
NSs and BHCs in quiescence (in the absence of other information). Section [| discusses the 
state of our observational understanding of the energy source for the quiescent emission from 
the NSs. We conclude in § [5] by summarizing and briefly discussing the application of this 
work to X-ray sources in globular clusters. 



2. Object Selection and Data Analysis 

Our purpose in fitting the H atmosphere model - appropriate only for NSs - to the 
BHC data is to directly compare the measured spectral parameters of the BHCs to those of 
NSs. We selected BHCs for analysis from among those in Table 1 of Menou et al. (|l9~9~9l) , 



which contains a list of 8 compact binary systems with implied masses ^3M . Three of 
these systems (A0620-00, GS 2023+33, and GRO J1655-40) were detected in X-rays in 
quiescence (at luminosities ^10 33 erg s^ 1 ), the data for which we analyse; the remaining 
five have upper limits in their luminosity. We analyse the three with luminosity upper-limits 
below 10 33 erg s" 1 (GRO J0422+32, GS 2000+25, and GS 1124-68) to investigate the 
constraints on emission area radius as a function of surface temperature. The remaining two 
sources with high luminosity upper-limits >10 33 erg s _1 (N Oph 1977 and 4U 1543—47) are 
consistent with or greater than the luminosities from NSs ( [Rutledge et al. 1999 ), and therefore 



they cannot be excluded as NSs based on a spectral comparison; we do not investigate them 
here. 

Brief descriptions of the analysed observations are in Table |I|. Data were obtained from 
the public archive at HEASARC/GSFC 0. All observations analysed are listed in Table [IJ 
all were performed with ROSA T/PSPC, except for those of GRO J1655— 40, which were 
performed with ASCA, and one observation of GRO J0422+32, which was performed with 
ROSAT/BBI. 

For our spectral fits and derived parameters, we adopted column densities using the 
conversion of iY" Hi 22 = A^ H /10 22 cm- 2 = 0.179 Ay flPredehl fc Schmitt l995|) and Ay = 



3.1 E(B — V) ( [Schild 1977| ), except where noted. We extracted the data for each source 
as described in the Appendix, and fit each extracted spectrum with a spectral model of 
galactic absorption and a tabulated H atmosphere model ( |Zavlin et al. 1996] ), using XSPEC 
( Arnaud 1996| ). 



The H atmosphere model is determined by two spectral parameters: an effective tem- 



8 



' http : / /heasarc . gsf c . nasa . gov/ 
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perature (kT eS ) and the ratio of an emission area radius (r e ; the true radius, which would be 
measured as the circumference of the NS, divided by 2n) to source distance. In our results, 
we quote r e by assuming a source distance (which can be uncertain by a factor of 2, and 
therefore represents a considerable systematic uncertainty in both the NSs and BHCs). The 
H atmosphere spectrum we used assumes the surface gravity of a 1.4M , 10 km object. It 
is possible to adopt the surface gravity as an additional parameter to the H atmosphere 
model; however, the available data are of insufficient quality to constrain simultaneously all 
3 parameters. 

We discuss in the Appendix the column densities adopted for each individual source. In 
quiescence, the S/N of the data is typically not sufficient to measure the X-ray equivalent 
column density for the assumed spectrum. We thus adopt an historically measured value 
of JVh, taken either from: (1) X-ray absorption observed while the object is X-ray bright; 
(2) optical reddening, which has been measured proportional to the equivalent hydrogen 
column density; or (3) the neutral hydrogen column density from radio observations ( Pickey 



fc Lockman 1990|) taken from the W3NH tool at HEASARC0, which measures the integrated 



column density not just to the distance of the X-ray object, but through the galaxy. All of 
these methods, when applied to estimating the JVh during X-ray quiescence, have systematic 
errors sufficient to produce >100% uncertainties in the implied emission area and surface 
temperature (see discussion in |Rutledge et al. 1999|) . In addition, there are observations 



suggesting a column density that varies over timescales of months to years, between outbursts 
of transients. An example is a change of AJVh = 0.5xl0 22 cm -2 in 4U 1608—522 QPenninx 



|et al. 1989D over several months. This uncertainty in A% can only be overcome by high 
S/N data in X-ray quiescence, which can permit a direct measurement of JVh during the 
observation. Interpretation of the results of the present analysis must bear this systematic 
uncertainty in mind. 



3. Comparison between BHCs and NSs in Quiescence 

In Fig. [| we compare the measured spectral parameters kT e g and r e of quiescent NSs 
and BHCs. Error bars are 90% confidence, as are upper-limits. For the four NSs in this 
analysis, the kT eS and r e were constrained, with best-fit kT e s in the range 0.08-0.20 keV, 
and r e in the range 8-12 km. The NSs exhibit a significant range in temperature, possibly 
due to different core temperatures (related to (M)). The data are of sufficient quality to 
constrain r e to within a factor of 2 for the NSs (not accounting for the uncertainties in 



'tittp : / /heasarc . gsf c .nasa.gov/docs/frames/mb_tools .html 
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source distance), and the resulting r e are consistent with objects of 10 km radius. The two 
connected points for Aql X— 1 are for two different distances (2.0 and 4.0 kpc). The two 
connected points for 4U 2129+47 are discussed in the analysis section [A.7| . We note that 



recent observations flUallanan ct al. 1999 ) have resolved the optical counterpart of Aql X— 1 



into two objects, only one (at most) will be associated with the X-ray system; most likely, 
this observation implies a distance to Aql X— 1 greater than previous estimates, although a 
new distance estimate has not yet been produced. 

The spectra of five of the six BHCs were of insufficient statistics to simultaneously 
constrain both r e and kT e s of the H atmosphere model. The exception (GS 2023+33) was 
significantly harder than the NS spectra, which constrained the kT e ^ to be above those 
observed from the NSs. The BHCs GRO J1655-40, GS 2000+25, GRO J0422+32 and 
GS 1124—68 have spectral parameters which are, within errors, consistent with those of the 
NSs. A0620— 00 and GRO J0422+32, if we presumed these to be NSs of the same surface 
area and gravity, must be cooler (^0.05 keV) than the average observed transient NSs. The 
spectrum of A0620— 00 is discrepant with those of the NSs, implying substantially smaller 
emission area (r e ~ 0.2-2 km) for the range of kT e g observed from the NSs. 

The addition of a power-law spectral component - such as the hard power-law tail 
as observed from Cen X-4 in ASCA data, and detected in Aql X-l in BeppoSAX data - 
introduces enough uncertainty in the H atmosphere spectral parameters of BHCs that all 
sources would be consistent with those we observe from the NSs. In noting this systematic 
uncertainty, we point out that additional spectral components are not demanded by the data 
we have used, although the data is largely of low-bandwidth (0.4-2.4 keV), and that wider 
bandwith data (such as from ASCA, Chandra, or XMM) may alter this. In practice, we 
have neglected this systematic uncertainty to investigate the question: if these BHC sources 
were discovered in quiescence, and their H atmosphere spectral parameters were compared 
with those from known NSs, would we conclude they are similar or dissimilar to the NSs? 
For GS 2023+33 and A0620-00, we find that they are dissimilar; while for GRO J1655-40, 
GS 1124-683, GS2000+25, and GRO J0422+32, we find they are consistent (albeit with 
a wide range of values as well). Higher signal-to-noise data, taken with an instrument of 
wider bandwidth, would permit us greater certainty regarding the possible contributions of 
spectral parameters for which we have not here accounted. 

Moreover, no BHC has well-constrained spectral parameters which would place it ex- 
clusively within the parameter space occupied by the NSs, as depicted on Fig. |2|. This is a 
phenomenologically defined region, which would have to be expanded with the discovery of 
NSs in quiescence outside of this box. Thus, we find no evidence that any of the BHCs have 
been misclassified, and should be re-classified as NSs. 
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4. The Energy Source for the Neutron Star Quiescent Luminosity 

There are presently no observational results which exclude that part of the quiescent 
luminosity of these NSs is due to accretion. Brown et al. (1998) noted a few observational 
tests which can be applied to determine whether accretion is active. First, as accretion 
will increase source luminosity (at this low luminosity level), the spectra should be drawn 
when sources are at their lowest observed luminosity in the X-ray passband. Second, since 
the H atmosphere thermal flux is expected to be variable only on timescales of longer than 
^months to 10 4 yr, variability on timescales less than this (<;days) likely indicates active 
accretion. Second, active accretion onto the NS surface will produce metal absorption lines in 
the spectrum (O and Fe, below 1 keV), which can be observed with Chandra and XMM, but 
which cannot be observed with ASCA or ROSAT. The presence of photospheric absorption 
metal lines in the spectrum - aside from being observationally important - will indicate 
active accretion onto the NS surface. These indicators should be used to insure that the 
observed NS emission is not due to active accretion (in addition to being observationally 
consistent with the theoretical H atmosphere spectrum). Since present instrumentation is 
not capable of detecting the lines, we can only apply the variability and observed luminosity 
criteria. 

If accretion is occuring during quiescence, this will increase the kT e g of an emergent H 
atmosphere spectrum, and produce metal absorption lines; accretion will not affect r e , unless 
emission originates from a different surface than the NS (such as in an accretion disk). 

Some observational evidence suggests that accretion is indeed occurring onto the NS 
surface during quiescence; long-term (months-years) variability in the observed flux has been 
reported (for 4U 2129+47, see App. |A.7.1| and |Garcia fc Callanan 1999| ; for Cen X-4, [Van 



[Faradijs et al. 1987|) . While this variability can be explained by a variable absorption column 



depth, active accretion during quiescence is also a possibility. However, recent observations 
of Aql X— 1 at the end of an outburst showed an abrupt fading into quiescence flUampana 



et al. 1998a|) associated with a sudden spectral hardening QZhang et al. 1998a]) . This was 



followed by a period of ~15 days, over which the source was observed (three times) with a 
constant flux level QCampana et al. 1998a|) . This behavior was interpreted as the onset of the 
"propellor effect" (|Illarionov fc Sunyaev 1975| ; |Stella et al. 1986|) in this object, which would 



inhibit - perhaps completely - accretion from the disk onto the NS. The energy source for 
the long-term nearly constant flux is then a puzzle. Unlike the accretion-powered models, 
the work of Brown et al. (1998) makes a specific prediction for this constant flux level, 
relating it to the long term time-averaged accretion rate. 



So as to minimize any contributions from accretion, we only analyse observations made 
during periods of the lowest observed flux. We calculate the bolometric luminosity from the H 
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atmosphere fits. The fits are given in terms of an (unredshifted) effective temperature (kT e g) 
and emitting area radius (r e ). The observed bolometric flux is F bol = (1 + z)~ 2 (r e /do) 2 aT^, 
where d is the source distance, and the observed bolometric luminosity is L bol = ind^F^. 
For a NS of 1.4 M and 10 km radius, the surface redshift is 1 + z = (1 - 2GM/Rc 2 )' l/2 = 
1.31. In Table § we give the unabsorbed, observed X-ray and bolometric luminosities as- 
suming a H atmosphere spectrum and the surface redshift of a 1.4 M 10 km radius NS 
(1 + z — 1.31), measured as described in the analysis of the present work and previously 
( IRutledge et al. 1999Q . The luminosity quoted for 4U 1608-522 in the previous work was 
not calculated for the distance stated there; the correct value is in this table. While there is 
considerable systematic uncertainty in these values, due to uncertain distances and A^h, the 
luminosities cover an order of magnitude, at about 10 32-33 erg s _1 . 

Using these new bolometric quiescent luminosities for Aql X-l, Cen X-4, and 4U 1608—522, 
we have remade a plot (Fig. |3|) of L q /L Q as a function of t r /t Q ( [Brown et al. 1998|) . Here 
L q and L Q are the observed quiescent and average outburst luminosities, and t r and t Q are 
the recurrence interval and outburst duration. We show this relation for the NSs (open 
circles) Aql X-l, Cen X-4, 4U 1608-522, and EXO 0748-676 and the BHCs (filled circles) 
H 1705-250, 4U 1543-47, Tra X-l, V 404 Cyg (GS 2023+33), GS 2000+25, and A 0620-00. 
We denote with an arrow those BHCs for which only an upper limit on L q is known. Be- 
cause the recurrence time for the NS 4U 2129+47 is unknown, we do not show it here. The 
expected incandescent luminosity is plotted for two different amounts of heat per accreted 
nucleon stored in the core during an outburst: lMeV(solid line) and 0.1 MeV(dotted line). 
The use of bolometric quiescent luminosities moved Aql X-l, Cen X-4, and 4U 1608—522 
to higher L q /L Q (upwards on this diagram). With the exception of these three objects and 
the Rapid Burster (L q is from [Asai et al. 1996a|) , the data from this plot is taken from 



( Chen et al. 1997| ). For Aql X-l and the Rapid Burster, L and t Q are accurately known 
(RXTE / All-Sky Monitor public data); for the remaining sources L Q and t D are estimated 
from the peak luminosities and the rise and decay timescales. 

Four of the five NSs are within the band where the quiescent luminosity is that ex- 
pected when the emitted heat is between 0.1-1.0 MeV per accreted baryon. The fifth NS 
(EXO 0748—676), has a higher quiescent luminosity (by a factor of 10), which we interpret 
as being due to continued accretion, an interpretation which is reinforced by the observation 
of spectral variability during the quiescent observations with ASCA ( |Corbet et al. 1994 ; 



Thomas et al. 1997| ), on timescales of ~1000 sec and longer. (Garcia and Callanan |1999 



measured L a .=lxl0 34 erg s _1 from Einstein/IPC observations of this source.) The BHCs 
on this figure are more spread out across the parameter-space, qualitatively indicating a 
statistical difference - although not one which is particular for each object - between the 
two classes of objects. This suggests the NS quiescent luminosity is more strongly related to 
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the accreted energy than the BH quiescent luminosities. 



5. Conclusions 

We have fit the quiescent X-ray spectra of four transient NSs and six transient BHCs 
(with measured mass functions) with a pure H atmosphere spectrum. We compared the 
emitting area radius and effective temperature of the ten sources and found that the NSs 
are clustered in (r e , kT e g) parameter space. Two of the BHCs (A0620— 00, GS 2023+33) are 
inconsistent with the NS spectra; the upper-limits of three more (GS 1124—68, GS 2000+25, 
GRO J0422+32), and the (r e , kT eS ) locus of the sixth (GRO J1655-40) overlap the NS pa- 
rameter space. The upper-limits of the parameter space of GS 2000+25 and GRO J0422+32 
are marginally consistent with the observed (r e , kT e s ) of the NSs, indicating that these BHCs, 
if interpreted as NSs, would have to be cooler than the average NS observed in quiescence. 

We found that the X-ray spectra excludes A0620— 00 and GS 2023+33 from being NSs 
of the type we used for comparison. These differences between the quiescent X-ray spectra 
of BHCs and NSs cannot be directly attributed to an observational selection effect. The 
NSs were identified by the type I X-ray bursts, the BHCs were first identified by their X- 
ray variability behavior and later by the high (^3M ) implied mass of the compact object. 
The X-ray spectra of GS 1124-68, GS 2000+25, GRO J1655-40, and GRO J0422+32 are 
consistent with those of the NSs, but their parameter space is large, and better data are 
needed to determine whether they are spectrally similar or not. 

Our method can be applied to objects which have low X-ray luminosities (<^10 34 erg 
s _1 ), to identify them as neutron stars, which have an atmosphere from which the theoretical 
emission originates, or as some object which is not a NS, as we have done here. In addition 
to the transient field objects, this method can be readily applied to the low-luminosity X- 
ray sources observed in globular clusters ( Hertz fc Grindlay 1983|) which are thought to be 
cataclysmic variables (|Cool et al. 1995| ; |Grindlay et al. 1995| ), but may also be transient 



neutron stars in quiescence ( [Verbunt et al. 1984|) . X-ray spectroscopic determination can 



identify these objects as NSs radiating thermal emission from the atmosphere, or imply a 
different origin for the emission. As discussed above and elsewhere ( [Brown et al. 199~8| ), the 
quiescent luminosities of these sources are set by the time average accretion rate. Thus, 
the low luminosity (10 31 erg s -1 ) X-ray sources in globular clusters, if they were transient 
neutron stars in quiescence, would have (M) « 2xl0 -13 M yr _1 ; comparing this to Aql X- 
1, with (M) « lxl0~ 10 Mq yr" 1 (estimated from the RXTE/ASM lightcurve history), and 
a mean outburst interval of ~ 200 days, the low-luminosity X-ray sources would have Aql- 
like outbursts with recurrence times of ~250 yr, assuming the quiescent luminosity is in 
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steady-state with accretion. The estimated number of such sources down to this luminosity 
level is highly uncertain, about 1-10 per globular cluster ( |Verbunt et al. 1995| ). For ~200 
GCs in the galaxy, one expects ~1-10 such transients per year. Even assuming only 1/4 of 
all GCs are within distance to detect an Aql X-l-like transient with the RXTE/ASM (for 
a peak RXTE/ASM detection countrate of 10 c/s; Aql X-l has an ASM peak countrate 
of ~30 c/s), this produces an expected transient discovery rate of 0.3-3 per year, which 
is consistent with or greater than the observed discovery rate with RXTE/ASM (of no 
sources over a four-year period). Stronger constraints could be made with more sensitive all- 
sky monitoring observations over longer time-baselines, or with a more tightly constrained 
luminosity function of the low-luminosity sources in globular clusters. 

Higher quality X-ray data from the coming X-ray spectroscopy missions (Chandra/ 'ACIS, 
XMM and ASTRO-E) will permit this analysis to be performed with greater accuracy, in 
particular by permitting the simultaneous measurement of the X-ray column density — a 
dominant systematic uncertainty. These will also provide the means to account for possible 
contributions due to a hard-power law component in the BHCs. For example, a 15 ksec obser- 
vation with Chandra of a source with photon power-law slope of 2, and luminosity of 2x 10 31 
erg s -1 , A"h,22=0.2 at lkpc would produce a spectrum which can be excluded as a H at- 
mosphere (or blackbody) with an (unconstrained) column depth, with probability=6xl0 -6 . 
Finally, the high spectral resolution and countrates of these instruments will permit a search 
for short timescale variability and photospheric metal absorption lines, which would indicate 
ongoing accretion during X-ray quiescence ([Brown et al. 19981) . 
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A. Sources 



The assumed distances (used to calculate luminosity) and adopted N H for the spectral 
fits performed here are listed in Table |2|; values for the NSs we analysed previously (Aql X— 1, 
Cen X— 4, 4U 1608—522) are in the previous reference flRutlcdge et al. 1999 ). The results 



of the spectral fits are presented in Table which contains: (1) the dataset number (cf. 
Table P; (2) the A^ H ; (3) the best-fit spectral parameters for the H atmosphere model, 
including the ■un-redshifted effective NS surface temperature (kT e g) and apparent emission 
area radius r e (km), as well as the reduced xt f° r that model. The H atmosphere model 
- appropriate only for NSs - is applied here to the BHCs, to produce spectral parameters 
which can then be directly compared with those observed from NSs. 

Here, we describe our assumptions about individual objects, and provide specifics of the 
analyses of each observation. 



A.l. GRO J0422+32 

The color excess, measured from an IUE spectrum of this object, was found to be 
E(B - V)=0.40±0.07 (]Vh,22~0.23 ±0.04; [Shrader et al. 1992|) . 

First, we analyse observation 1, a ROSAT/PSPC observation which had not previously 
been analysed (this observation is included for completeness, as this source was fainter during 
observation 2). We extracted the data from a 53" circle about the source, and background 
from a 530" annulus about the source, excluding another object in the FOV, and 3 low-surface 
brightness areas clustered in the SE side of the annulus (possibly background fluctuations). 
Due to the long exposure, this spectrum constrained all three parameters (Ah, /cT eff , and 
r e ). The resulting best-fit was r e =3.0tj$ km, fcT eff =0.22lg;^ keV, and A H , 2 2< 0.12 (90%), 
for an implied luminosity of 7.8xl0 32 (d/2.6 kpc) 2 erg s _1 (0.5-2.0 keV). In a fit holding the 
column density fixed at ^,22=0.22, the best fit model was not acceptable at the p=1.8% 
level. The parameter space for a A% 2 =2.72 would constrain the r e to be similar to those 
found from the NSs (llif^ km). 

We analysed the HRI observation 2. The observation was previously analysed, pro- 
ducing a luminosity upper-limit of log(L)<31.6 (d=2.6 kpc, a=2.1, and Ah,22=0.2; |Garcia 
|et al. 19971) . We extracted the spectrum from an 8" radius about the source position - finding 
7 photons, consistent with the expected and measured background count-rate. We extracted 
background counts from an annulus of 50" and 10" outer and inner radii, respectively. We 
used the Dec 1 1990 HRI spectral response from GSFC/HEASARC calibration database. 
We find a slightly lower, but consistent 3a upper-limit to the unabsorbed source luminosity 
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for the same assumed spectrum (log(L) < 31.4) as found previously. We assumed a series 
of temperatures and found upper-limits to r e for the assumed spectrum. For a comparable 
temperature as derived from the PSPC analysis (above), the implied r e < 0.8 km - sub- 
stantially below that found from the PSPC analysis, which indicates intensity variability 
at this low luminosity. As observation 2 has lower luminosity than observation 1, we use 
observation 2 in our interpretations (as thermal emission from the NS surface should not 
vary by ^few% over the timescales here, and it is only the lowest observed luminosity from 
each object which may be due in its greatest part to the thermal surface emission). 



A.2. A0620-00 

The optical color excess of the companion star has been measured as E(B—V) =0.39 ±0.02. 
( |Wu et al. 1976| ; pke fc Greenstein 1977| ), corresponding to A" H ,22=0.22 - somewhat below 
the ./Vh, 22=0. 41 (pickey fc Lockman 1990| ) measured in this direction (which is the integrated 



value along this line of sight through the entire galaxy), consistent with a nearby object. We 
adopt the ^,22=0.22 value for our spectral fits. 

We analysed data taken by ROSAT /PSPC (observation 3); this data was previously 
analysed ( |McClintock et al. 19951) . An X -ray source was found in the extracted data at 
the source position. Counts from this source were extracted from a circle of radius 53". 
Background was taken from an annulus about the source, 4' outer radius, and 60" inner 
radius; excluded from this annulus were circular areas (each only half-within the annulus, 
straddling the outer radius) about two unrelated sources, each with radius 53". The source 
region contained a total of 116 counts. The background region contained 938 counts. We used 
data in the energy range 0.4-2.1 keV, rebinned into three energy bins 0.4-0.8, 0.8-1.2, and 1.2- 
2.1 keV. Using the same spectral model and assumptions as the previous work ( |McClintock 



et al. 199"5| ), we reproduce the best-fit blackbody temperature and source luminosity. 



The data are of insufficient S/N to constrain the temperature and area independently. 
We held the temperature fixed at values kT e s= [0.04, 0.54] keV, and extracted the best 
fit emission area radii (r e ), which were found to range from 14 ±2 km (at kT c ^= 0.04 
keV) to 0.07 km (at kT eS = 0.15 keV). The model xl reaches a probability of 2% that the 
observed spectrum is produced by the model in a single random trial (xt = for 2 dof) 
at kT e {[= 0.265 keV, which sets our upper-limit on the effective temperature for this source. 
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A.3. GS 1124-68 

The X-ray measured equivalent column density was found to have evolved, decreasing by 
-^h,22~0.2 during the outburst, and settling to a value of A^ H ,22=0.16 ( [Ebisawa et al. 19"§4| ), 
consistent with the optical reddening E{B — y)=0.20±0.05 flGonzalcz-Riestra et al. 1991 ) 



measured from a UV absorption line, and is below the measured hydrogen column density in 
the direction of this source ./Vh, 22=0. 25 flDickey fc Lockman 1990|) , consistent with an object 



a short distance away relative to the size of the galaxy. We adopt ./Vh,22=0.16 for our spectral 
fits. 



The PSPC observation (number 4) has been previously analysed ( |Greiner et al. 1994 



INarayan et al. 1997b[). For the source, we extracted counts from a 75" radius circle about the 



source position, as the source position is about 15' off-axis. For background, we used three 
circular areas on the detector, one centered at the source position, the other two centered 
at distances from the FOV center equal to that of the source position. From these three 
circular areas, we excluded the 75" region about the source position, and one other apparent 
(that is, faint) source. There were 66 counts in the source region, and 2104 counts in the 
background spectrum. We generated an ancillary response file for this off-axis spectrum 
using the FTOOL pcarf V2.1.0. In the spectral analyses, we find a 3a upper-limit to the 
(unabsorbed) source flux of <3.5xl0~ 14 erg cm -2 s -1 (0.3-2.4 keV), assuming iV*H,22 = 0.22, 
and a power-law photon spectral slope of a = 2.5, consistent with the previously found 
value. 

There is no detectable quiescent source at the position of the optical source, consistent 
with previous results ( |Greiner et al. 1994| ; [Narayan et al. 1997b| ). Across the investigated 



range of NS surface temperatures, the upper-limits of the implied radii range from 82-0.13 
km. 



A.4. GRO J1655-40 

This ASCA observation has been analysed previously ( |Ueda et al. 1998] ), in which the 
flux was measured to be 2xl0~ 13 ergcm~ 2 s _1 (2-10 keV). The interstellar column density 
was measured at A^h,22 = 0.74±0.03 in the X-ray high (bright) state, and <0.08 and <0.14 
(90%)in the low (faint) state. The optical reddening is E(B — y)=1.3±0.1 (cf. Prosz fc] 
Pailyn 1997| ; |Horne et al. 1996|) is consistent with a value of ./Vh,22=0.74, which we adopt 



and hold fixed. 



Using data taken with the GIS in PH mode, we extracted background counts from two 
circles, each 5' in radius and centered 5' from the image center - the same distance as the 
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position of the object. The source X-ray spectra were extracted from a 5' circle about the 
optical position. We excluded energy bins below 1.0 keV and above 5.0 keV from the fit. 

For the SIS, data were taken from 2.5' radius circles about the source. Background was 
taken from a square approximately 6' on a side, excluding the 2.5' source region. We used 
energy range of 1.0-10.0 keV for the GIS data, and 0.5-10.0 keV for the SIS data. 

Based on this data, no single best fit is found for 0.04keV <kT c g< 0.54keV. As the 
assumed temperature (held fixed) is increased from fcT eff =0.04 to 0.54 keV, r e decreases from 
110±32 km to 0.19±0.02 km. 



A.5. GS 2000+25 

The X-ray equivalent Ah for this object has been measured as Ah,22 = 1-14 Qlsunemi 
|et al. 1989]) , and optical reddening measured as E(B — V) = 1.5 ( phevalier fc llovaisky 1990]) 
(corresponding to Ah,22 = 0.82 5), which are slightly discrepant from one another, though 
within the systematic uncertainties between these techniques. The measured hydrogen col- 
umn density in the direction of this object is Ah,22=0.66 ( Pickey fc Lockman 1990|) , although 
within a degree, values range between A Hj 22= [0.51,0.81], and at the most proximate point, 
Ah,22=0.67. The higher measured optical reddening of the stellar companion and absorption 
of X-rays from the compact object, relative to the direct Ah measurement, may indicate 
significant absorption in the local environment of the binary, on the order of Ah i 22=0.5. We 
use in turn A Hj 22=0.66 and A Hi 22=l-l. 

Two of the presently analysed observations (#6 & 7) have been analysed previously 
( |Verbunt et al. 1994]) . In the ROSAT image for any of the three observations (#6, 7, and 8) 
there is no detected object at the source position. We combine the data, which are of three 
widely varying epochs, to produce average upper limits to spectral parameters. 

In each of the three observations, we extract source counts from a circular region centered 
on the object of 60" radius. Background counts were taken from annuli about the source, 
with an inner radius of 63.75" and an outer radius of 375". We also excluded two circular 
regions (60" radius) which overlapped the annuli, which appear to contain sources. There 
were a total of 78 counts in the source region, and 2505 counts in the background region. The 
resultant spectrum produces flux upper-limits consistent with (but below, by about 25%) 
the higher upper-limits found previously ([Verbunt et al. 1994j) . 

We do not detect a source above the background countrate; we produce 90% upper 
limits to r e , which were found to range from 0.13-28 km, depending on the assumed surface 



temperature, in the range of 0.04-0.54 keV. 



A.6. GS 2023+33 



The optical reddening has been measured at Ay ~ 3.0 ( Charles et al. 1989 ), and a 



comparable E(B - V)=1.03 flWagner et al. 199 ID , which imply A% >22 =0.54-0.59. This is 



below the measured A%=0.81 in this direction ( pickey fc Lockman 1990| ) consistent with an 
object nearby relative to the size of the galaxy. We adopt a value of A^h,22=0.54. 

This object was observed twice in quiescence: in 1994 with ASCA and in 1992 with 



ROSAT /PSPC (see |Verbunt et al. 1994]) . It was also observed and analysed in the ROSAT 
All-Sky Survey ([Verbunt et al. 1994|) , on Nov 5, with an average countrate of 0.028±0.009 
c/s,in which it was found with a luminosity of about 1.5xl0 33 (d/3.5 kpc) 2 erg s -1 (0.4-2.4 
keV). 

From the ROSAT image, data were selected centered on the X-ray source, in a circle 
60" in radius. Background was taken from an annulus with an inner-radius of 63.8" and an 
outer radius of 375". The source region had 425 counts, and the background region had a 
total of 1365 counts. 

At effective temperatures below, kT e g= 0.258, the model becomes untenable (x 2 ,=4.0, 2 
dof, corresponding to 98% probability), which sets the lower-limit of the acceptable effective 
temperature parameter space for this data. The best fit area was not well constrained from 
below (as the temperature increases to values above those for which this model is valid), 
but the 90% upper limit on the radius is 1.9 km (d/3.5 kpc). We find the (0.4-2.4 keV) 
luminosity of the best fit model to be 6xl0 32 (<i/3.5kpc) 2 erg s _1 . 



These results are consistent with those of an earlier analysis of these same data ([Wagner 



et al. 199"4| ), although we note that intensity variability was found on ~days timescales which, 



a priori cannot be attributed to the NS surface flux from the hot (thermal) core. These results 
are also consistent with results from ASCA analysis of quiescent data taken at a different 
time ( Narayan et al. 1997a| ). 



A.7. 4U 2129+47 



We inadvertently did not include this NS transient in our previous study, so we analyse 
it here. This spectrum has been analysed previously with a black-body model flGarcia fc| 
Callanan 1999|) , and the derived parameters were consistent with those we find here, under 
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similar distance /column density assumptions. 

The distance to this object is controversial with values between 1 and 6 kpc (see |(Jowley 



fc Schmidtke 1990| for arguments on both sides). Associated with this distance uncertainty 
is an optical reddening uncertainty; if the optical counterpart is more distant, it must be 
more luminous, of earlier type, and bluer, in which case the reddening is E(B — V)=0.3 
(7X^22=0. 17; |Cowley fc Schmidtke 1990| ); if less distant, the companion is redder, and EiB — 
V)=0.5 (A^ H ,22 =0.28; |Thorstensen et al. 1979| ; |Ghevalier et al. 1989Q . The total N H measured 



in the direction of this object is A r n,22 =0.38 (pickey fc Lockman 1990|) , while the best fit 



column density obtained by Garcia and Callanan (with this same data; |1999| ) was Ah,22=21.5 
(no quoted error). 

In interpreting the data, we alternately adopt a distance of 1.5 kpc and EiB — V)=0.5 
(AT H 22 =0.28), and 6.0 kpc and E(B - V)=0.3 (N R , 22 =0.17). 

For the March 1994 PSPC observation, the source counts were extracted from a circular 
region 45" in radius, with a total 279 counts. Background was extracted from an annulus 
about the source, 278" and 45" in outer and inner radius, respectively, for a total of 1959 
counts. We fit to data in the 0.4-2.4 keV energy range. 

For an assumed ]Vh,22=0.28 (d=1.5 kpc), the best fit parameters were kT cS =0.08tox)i 5 
keV, and r 6 =7.1±|;g km. For jV Hj2 2=0.17 (d=6.0 keV), the best fit parameters are kT eS =0.l05±° ° f , 
and r e =12.0^4;g. The unabsorbed luminosities for the assumed (A%, d) combinations are 
6.5xl0 31 and 5.6 xlO 32 erg s _1 (0.4-2.4 keV), respectively. 



A. 7.1. An Observed Change in the Flux from 4U 2129+47 

To compare the observed flux of 4U 2129+47 analysed here, with one measured pre- 
viously with the ROSAT /HRI during a 7 ksec of observations (Nov-Dec 1992; |Garcia 



|et al. 199TD , we imposed a A^"h,22=0.5 and a photon power-law slope of 2.0, and mea- 
sured the (absorbed) flux from the March 1994 observation, which was (6.7±0.5) x 10~ 14 
erg cm -2 s _1 (0.3-2.4 keV) - a factor of 3.4±0.6 below that measured during the earlier ob- 
servation by Garcia (we independently analysed the earlier HRI observation, and confirm 
Garcia's flux measurement). We cannot simultaneously extract r e and kT e g from the HRI 
observation, as only one spectral bin is available, and so we cannot compare the photon en- 
ergy spectra. However, we jointly fit both observed spectra with the same H atmosphere and 
A^h, and found that there is a 2.7% chance that the best-fit spectrum produced both observed 
spectra. Assuming that r e and A^h are the same between the two observations, the effective 
temperature would have to have dropped from 0.091±0.006 keV (90%) to 0.081±0.002 keV 
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to explain the different spectra. The thermal timescale at the depth where the crust heating 
occurs is about a year ( [Brown et al. 1998[ ); the decrease in effective temperature could be a 
signature of a cooling crust. Another possible explanation is that the Nn was lower during 
the 1992 observation than the 1994 observation; for example if A^h,22=0.10 during the HRI 
observation and ]Vh.22=0.28 during the 1994 observation, but with the same underlying H 
atmosphere spectrum, this would account for the observed difference in flux. In addition, if 
accretion is active during quiescence, the different X-ray spectra may reflect a decrease in 
the accretion rate by a factor (0.081/0.0906) 4 =0.64, assuming a thermal spectrum. 

With the quality of spectra we have here, we cannot definitely state the cause of the 
discrepancy between the 1992 and 1994 observations, but it is consistent with either a differ- 
ence in the intervening A%, a difference in the effective temperature, time- variable accretion, 
or a combination of the three. 
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Fig. 1. — Comparison between the spectral parameters r e and kT e g, derived from spectral fits 
of the quiescent X-ray emission from Aql X-l, Cen X-4, 4U 1608-522 and 4U 2129+47. 
The (open points) are for the H atmosphere spectrum and the (solid points) are from a 
black-body spectrum. The two points connected for Aql X— 1 correspond to the upper- and 
lower- distance limits for that source. The two connected points for 4U 2129+47 are for the 
two different distance/A^ estimates (see text). As expected from comparisons between the 
theoretical spectra, the H atmosphere fits produce higher values of r e . 

Fig. 2.— The best fit emission area radii for NSs and BHCs (A0620-00, GS 2023+33, 
GRO J 1655— 40) using a hydrogen atmosphere spectrum, with the surface gravity of a 1.4M , 
10km object. The dotted-line box encloses the (phenomenological) parameter space occupied 
by the NSs. The 90% upper-limits for GS 2000+25, GS 1124-68 and GRO J0422+32, are not 
included (see Table |D, but are lines which lie between GRO J1655— 40and A0620— OOacross 
kT eS =0. 04-0. 054 keV, passing through the NS box. The best-fit parameters for the NSs 
cluster in the /cT e g~0.08-0.2 keV and r e ~8-12 km. The best-fit parameters for BHCs 
GRO J1655-40, GS 1124-68, GS 2000+25, and GRO J0422+32 are consistent with having 
similar spectra, while GS 2023+33 and A0620— 00 are inconsistent. 
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Fig. 3. — The ratio of outburst luminosity L q to quiescent luminosity L a as a function of 
the ratio of recurrence interval t r to outburst duration t . The lines are for different amounts 
of heat, 0.1 MeV (dashed line) and 1.0 MeV (solid line), per accreted nucleon deposited at 
depths where the thermal time is longer than the outburst recurrence time. Also plotted 
are the observed ratios for several NSs (filled squares) and BHs (open squares). For most 
of the BHs, only an upper limit (arrow) to L q is known. Data is from ( |Chen et al. 1997| ), 
with the exception of L q for the Rapid Burster ( [Asai et al. 1996151) . For Aql X-l and the 
Rapid Burster, L and t a are accurately known (RXTE / All-Sky Monitor public data); for the 
remaining sources L and t Q are estimated from the peak luminosities and the rise and decay 
timescales. For Aql X-l, Cen X-4, and 4U 1608—522, we have recalculated the bolometric 
luminosities using the H atmosphere models ( |Rutlcdge et al. 1999Q . 
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Table 1: Observations 



Number 


Satellite/Instrument 


Obs Start Time 


Live Time 


Avg Countrate a 






(UT) 


(ksec) 


(c/s) 






GRO J0422+32 






1 


ROSAT/PSPC 


17 Aug 1993 04:37 


28.0 


(8.0+0.2) xl0~ 2 


2 


ROSAT/ma 


6 Jun 1995 12:29 


18.8 


(3.3+3. 3)xl0~ 4 






A0620-00 






3 


ROSAT /PSPC 


10 Mar 1992 17:05 


29.8 


(1.5±0.3)xl0~ 3 






GS 1124-68 






4 


ROSAT /PSPC 


1 Mar 1992 00:11 


16.0 


(4.0+3. 5)xl0~ 4 






GRO J1655-40 






5 


ASCA 


23 Mar 1996 00:54 


87.0 


(2.1+0.4) xl0~ 3 (SIS1) 






GS 2000+25 






6 


/20&4T/PSPC 


1 May 1992 13:23 


12.4 


(0.8±6.0)xl0~ 4b 


7 


EO&4T/PSPC 


9 Apr 1993 06:41 


6.6 




8 


ROSAT/PSPC 


13 Oct 1993 13:11 


6.8 








GS 2023+33 






9 


ROSAT 


5 Nov 1992 13:04 


16.5 


(2.1±0.1)xl0~ 2 






4U 2129+47 






10 


ROSAT /HRI 


9 Dec 1991 12:47 


7.0 


(4.8±1.0)xl0- 3 


11 


EO&4T/PSPC 


6 Mar 1994 00:10 


31.0 


(6.7+0.5) xl0~ 3 



a Count-rates are background subtracted. For ROSAT/PSPC and HRI, countrates are for the energy range 
0.4-2.4 keV; for ASCA the relevant instrument is listed. 
^Average for observations 6, 7, & 8 



Table 2: Object Information 



Object 


do 


Mass 


N u 




(kpc) 


(M e ) 


(10 22 cm" 2 ) 


Black Hole Candidates with Mass Estimates 


GRO J0422+32 


2.6 1 


>9 2 


0.22 3 


A0620-00 


0.9 4 


>7.3 5 


0.22 6 ' 7 


GS 1124-68 


3.3 8 


5.0-7.5 9 


0.16 10 


GRO J1655-40 


3.2 11 


7.02+0.22 12 


0.70 13 < 14 < 15 


GS 2000+25 


2.0 16 


7.0-7.7+0.5 17 


0.66-1.14 13 ' 16 ' 18 


GS 2023+33 


3.5 19 


8-12 19 


0.57 20 




Neutron Star 




4U 2129+47 a 


1.5 21 


(NS) 


0.28 21 


4U 2129+47 6 


6.0 22 


(NS) 


0.17 22 



References. — l,[Esinctal. 199S| ; 2, [Beekman et al. 1997|; 3, |5hrader et al. 1992| ; 4, |Okc fe Grccnstcin 197^ ; 
5, |McClintock fe Remillard 1986|; 6, |Wu et al. 1976|; 7, pke fc Greenstein 1977|; 8, |Shahbaz et al. 1997] ; 9, 



Orosz ct al. 1996|; 10, |Ebisawa et al. 1994j 11, |Hjcllming fc Rupcn 1995 ; 12, |Orosz fc Bailyn 1997|; 13, W3 nH, 
Dickey fc Lockman 1990| ; 14, |Ueda et al. 1998j 15, |Horne et al. 199§ 16, |Chevalier fc Ilovaisky 1990| ; 17, 



et al. 1979 



Casares et al. 1995 ; 18, rsunemi et al. 1989] ; 19, Wagner et al. 1992 ; 20, |Wagncr et al. 1991 ; 21, [Thorstensen 
22, powley fc Schmidtke 1990| ; a, |Thorstensen et al. 1979) phcvalicr et al. 1989| ; b, powlcy fc 



Schmidtke 1990 

Gonversion factors: 7V H ,22 =0.179 A v ; ,4^=3.1 E(B - V). See Text 



Table 3: Derived Spectral Parameters 



Datasct 




fcT off 




X^(dof) 


Number 


(10 22 cm -2 ) 


(kcV) 


(km (d/d )) 








GRO .10422+32 






1 


<0.12 


0.14+0.01 


, + 2.4 


1.64/1 




(0.22 } 


U. 14±U.U1 


1 1 + 4 - 7 
n -2.5 


6.3/2 


2 


(0.22) 


(0.04) 


<27 


n/a 




(0.22) 


(0.08) 


<3.6 


n / a. 




(0.22) 


(0.15) 


<0.8 


n / a. 




(0.22) 


(0.30) 


<0.2 


n / a. 




(0.22) 


(0.54) 


<0.08 


n/a 






A0620-00 






3 


(0.22) 


(0.04) 


14+2 


1.2/2 




(0.22) 


(0.08) 


1.8 ±0.3 


0.3/2 




(0.22) 


(0.15) 


0.37 ±0.06 


1.9/2 




(0.22) 


(0.265) 


0.105 


4.0/2 






GS 1124-68 






4 


(0.16) 


(0.04) 


<82 


n / a 




(0.16) 
(0.16) 


(0.08) 
(0.15) 


<8.6 
<1.2 


ii /a 

n /a 




(0.16) 


(0.30) 


<0.26 


n/a 




(0.16) 


(0.54) 


<0.13 


n/a 






GRO J1655-40 






5 


(0.74) 


(0.04) 


110±32 


105/70 




(0.74) 


(0.08) 


12.0±2.8 


100/70 




(0.74) 


(0.15) 


2.35±0.4 


95/70 




(0.74) 


(0.30) 


0.53±0.08 


80/70 




(0.74) 


(0.54) 


0.19±0.02 


61/70 






GS 2000+25 






6, 7, & 8 


(0.66) 


(0.04) 


<28 


n/a 




(0.66) 


(0.08) 


<3.8 


n/a 




(0.66) 


(0.15) 


<1.0 


n/a 




(1.1) 
(0.66) 


(0.15) 
(0.30) 


<1.9 

<0.29 


n/a 
n/a 




(0.66) 


(0.54) 


<0.13 


n/a 






GS 2023+33 






9 


(0.54) 
(0.54) 


(0.258) 
(0.30) 


2.7+0.1 
2.0+0.1 


4.0/2 
1.8/2 




(0.54) 


(0.54) 


0.81±0.03 


0.72/2 


11 


4U 2129+47 (1.5 kpc; AT H 
(0.28) 0.08t° ° 2 15 


22=0.28) 

71 +6.0 
— 3 5 


0.13/1 


11 


4U 2129+47 (6.0 kpc; JV H 
(0.17) 0.105+2;^ 


22=0.17) ' 
19+8.6 
iZ -4.9 


0.71/1 



Errors and upper-limits are 90% confidence; values listed in parenthesis were fixed during the fits. For 
description of the different (do, Nn) assumptions for 4U 2129+47, see App. A. 7 



Table 4: Bolometric Luminosity of Transient Neutron Stars in Quiescence 



Source L x (0.5-10.0 keV) L bo i 

10 32 erg s _1 10 32 erg s~ 



Aql X-l 




(20.4+9.2) (V4kpc) 2 


39 


Cen X-4 




(1.6+0.6) (rf /l-2kpc) 2 


4.6 


4U 1608-522 




(30.0+15.2) (rf /3.6kpc) 2 


56 


4U 2129+47 a 


(0.7 H _ 


:H) (rf /1.5kpc) 2 (0.5-3.0 keV) 


1.6 


4U 2129+47 b 


(5.6; 


n 3 2 ) (do/6.0 kpc) 2 (0.5-3.0 keV) 


13 



■ b A range of (ATjj, do) is assumed; See App. 



A.7 



